Study of Liapunov Exponents and the Reversibility of Molecular Dynamics Algorithms

نویسندگان

  • Karl Jansen
  • Chuan Liu
چکیده

Molecular dynamics algorithms like the Hybrid Monte Carlo [1] and the Kramers equation [2,3] algorithm have been playing a major role in numerical simulations of QCD on a Euclidean lattice. In order for these algorithms to fulfill the detailed balance condition, the classical motion governed by the set of Hamilton’s equations for the system should be reversible. However, as noticed some time ago [4], this reversibility condition is violated due to the round-off errors in the numerical integration of the equations of motion. Recently, it was pointed out [3] that, due to the chaotic nature of the equations of motion, the round-off errors in the integration which violate the reversibility condition get magnified exponentially with a positive Liapunov exponent ν. We will focus our discussion on the equations that arise in the simulations of Wilson QCD with two flavors of quarks with degenerate masses. Periodic boundary conditions have been taken for all the fields. The full partition function for Wilson QCD is given by,

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Instabilities and Non-Reversibility of Molecular Dynamics Trajectories

The theoretical justification of the Hybrid Monte Carlo algorithm depends upon the molecular dynamics trajectories within it being exactly reversible. If computations were carried out with exact arithmetic then it would be easy to ensure such reversibility, but the use of approximate floating point arithmetic inevitably introduces violations of reversibility. In the absence of evidence to the c...

متن کامل

Energy study at different solvents for potassium Channel Protein by Monte Carlo, Molecular and Langevin Dynamics Simulations

Potassium Channels allow potassium flux and are essential for the generation of electric current acrossexcitable membranes. Potassium Channels are also the targets of various intracellular controlmechanisms; such that the suboptimal regulation of channel function might be related to pathologicalconditions. Realistic studies of ion current in biologic channels present a major challenge for compu...

متن کامل

Liapunov Spectra for Infinite Chains of Nonlinear Oscillators

We argue that the spectrum of Liapunov exponents for long chains of nonlinear oscillators, at large energy per mode, may be well approximated by the Liapunov exponents of products of independent random matrices. If, in addition, statistical mechanics applies to the system, the elements of these random matrices have a distribution which may be calculated from the potential and the energy alone. ...

متن کامل

Asymptotic Liapunov Exponents Spectrum For An Extended Chaotic Coupled Map Lattice

The scaling hypothesis for the coupled chaotic map lattices (CML) is formulated. Scaling properties of the CML in the regime of extensive chaos observed numerically before is justified analytically. The asymptotic Liapunov exponents spectrum for coupled piece-wise linear chaotic maps is found.

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 1996